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Density-functional theory study of the initial oxygen incorporation in Pd(lll) 
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Pd(lll) has recently been sliown to exliibit a propensity to form a sub-nanometer tliin surface 
oxide film already well before a full monolayer coverage of adsorbed O atoms is reached on the surface. 
Aiming at an atomic-scale understanding of this finding, we study the initial oxygen incorporation 
into the Pd(lll) surface using density-functional theory. We find that oxygen incorporation into 
the sub-surface region starts at essentially the same coverage as formation of the surface oxide. This 
implies that the role of sub-surface oxygen should be considered as that of a metastable precursor 
in the oxidation process of the surface. The mechanisms found to play a role towards the ensuing 
stabilization of an ordered oxidic structure with a mixed on-surface/sub-surface site occupation 
follow a clear trend over the late 4d transition metal series, as seen by comparing our data to 
previously published studies concerned with oxide formation at the basal surface of Ru, Rh and Ag. 
The formation of a linearly aligned O-TM-0 trilayered structure (TM = Ru, Rh, Pd, Ag), together 
with an efficient coupling to the underlying substrate seem to be key ingredients in this respect. 

PACS numbers: 81.65.Mq, 68.43.Bc, 68.47.Gh 
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I. INTRODUCTION 

Contact with a rich gas-phase is characteristic for solid 
surfaces in most applications and real-life situations. The 
concomitant interaction with the atoms and molecules of 
the environment can induce significant materials changes, 
exemplified prominently by the corrosion of metal sur- 
faces. The oxide formation underlying the latter is an 
everyday occurrence in our oxygen-rich atmosphere and 
of great importance for the atomic-scale understanding 
of many technologically relevant processes like oxidation 
catalysis. There is general agreement that the atomistic 
steps leading to oxide formation involve dissociative oxy- 
gen chemisorption on the metal surfaces, often the forma- 
tion of adsorbate structures, followed by incorporation of 
oxygen into the surface and finally (or simultaneously) 
the transformation to oxidic films (surface oxides) and 
bulk oxide growth. Still, the microscopic understand- 
ing of these different processes is quite shallow. The se- 
quence of atomistic events may be quite complex, and 
rather than successively, they may occur simultaneously 
or even concertedly, depending crucially on temperature 
and partial pressures in the surrounding. 

Investigations on well-defined single-crystals help to 
improve this atomic-scale insight, and due to its 
widespread technological use, Pd is an appropriate tar- 
get material. The initial chemisorption of oxygen on 
Pd single crystal surfaces has correspondingly been the 
subject of i ntensive experimental research, e.g. Refs. 
I il l 2I3I I4II 51 16171 1 81. The picture arising from these studies for 
the close packed (111) surface is that for temperatures 
T > 200 K oxygen adsorption on the surface is dissocia- 
tive. An ordered (2 x 2) overlayer is observed at a cover- 
age of = 0.25 monolayers (ML)^, while at higher cov- 
erages on-surface adsorption seems to compete with the 
formation of oxidic structures: depending on the prepa- 
ration conditions a (\/3 x ^3) i?30°\ a (1 x l)^-^ and 



a "complex" structure^'* have been observed, where the 
"complex" structure with Q — 0.7 ML was recently char- 
acterized as a pseudo-commensurate surface oxide con- 
sisting of an 0-Pd-O trilayered geometry on the Pd(lll) 
substrate^. 



In a previous publication^ we discussed the on-surface 
adsorption of oxygen on the Pd(lll) surface for cov- 
erages up to one monolayer. Using density-functional 
theory (DFT) we identified the fee hollow sites as the 
most stable adsorption sitea^, which is consistent with 
a trend observed for the basal surface of the other late 
Ad transition metals, namely Ru, Rh and Ag, for which 
oxygen also adsorbs preferentially in those sites follow- 
ing the bulk-stacking sequence. Furthermore, we found 
that overall repulsive lateral interactions in the addressed 
coverage range give rise to a pronounced reduction in the 
adsorption energy with increasing oxygen coverage, but 
are still clearly exothermic up to a full ML (1 x 1) phase. 
In the present paper we now investigate the initial oxy- 
gen incorporation into the palladium (111) surface us- 
ing DFT. We find that occupation of sub-surface sites 
becomes energetically more favorable than continued ad- 
sorption into on-surface sites already at total coverages 
well below IML. We analyze the ensuing structures and 
discuss what might be the mechanisms which make one 
structure more stable than another. A comparison to the 
4d neighbors of Pd in the periodic table shows that the 
mechanisms found here are valid there, as well. Last, but 
not least, we address the question of how much oxygen 
can be dissolved into a palladium crystal, an issue which 
has often been the object of speculation in the experi- 
mental literature. 
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II. THEORY 

The DFT calculations are performed within the full- 
potential linear augmented plane wave (FP-LAPW) 
methodiSiiiii^ using the generalized gradient approxi- 
mation (GGA)i^ for the exchange-correlation functional. 
All surface structures are modeled in a supercell geome- 
try, employing a symmetric slab consisting of seven (111) 
Pd layers. A vacuum region of 17A ensures the decou- 
pling of consecutive slabs. All atomic positions within 
the adsorbed oxygen layers and the two outermost Pd 
layers are fully relaxed. 

The FP-LAPW basis-set parameters are as follows: 
the mufhn-tin radii are i?^T ~ 2.25 bohr and -R^t ~ 
1.3bohr. A wave function expansion inside the muffin 
tin spheres up to ZJJ^^x = 12 and a potential expansion up 
to = 4 is used. The 4s and Ap semicore states of Pd, 
as well as the O 2s states are described by adding local 
orbitals to the FP-LAPW basis set. The energy cutoff 
for the plane wave representation in the interstitial re- 
gion between the muffin tin spheres is E^^-^ = 17Ry for 
the wavefunctions and Ep]^{^ = 169 Ry for the potential. 
A (12x12x1) Monkhorst-Pack grid with 19 k points in 
the irreducible wedge is employed for the Brillouin zone 
integration in a (1 x 1) cell. To obtain the same sam- 
pling of the reciprocal space for bigger surface cells, this 
number is reduced accordingly. 

This is the same basis set as in our preceding publi- 
cation on the oxygen on-surface adsorption on Pd(lll)S. 
Its high accuracy has already been described there, and 
with respect to the present work it is primarily charac- 
terized by a ±30meV (±50meV) numerical uncertainty, 
when comparing relative binding energies in geometries 
containing an equal (unequal) amount of O atoms. 

III. RESULTS 
A. Initial oxygen incorporation 

This subsection addresses the initial oxygen incorpora- 
tion into the surface by studying the stability and prop- 
erties of a set of structures with oxygen coverages in the 
range < 9tot < IML. We employ (2 x 2) unit cells and 
calculate the average binding energy per O atom with 
respect to the dissociation energy of an oxygen molecule. 
Our sign convention is such that positive binding ener- 
gies indicate that adsorption is exothermic. From one 
(Btot = 0.25 ML) up to four (Gtot = 1.00 ML) oxygen 
atoms are placed in any combinatorial way that is possi- 
ble in a (2 X 2) unit cell in on-surface sites and in inter- 
stitial sites between the first and second substrate layers. 
Occupation of sub-surface sites will then obviously com- 
mence when a structure containing oxygen in sub-surface 
sites becomes energetically more stable than its counter- 
part with the same total adsorbate coverage, 9tot, but 
all oxygen atoms located on the surface. 

In our previous publication on oxygen on Pd(lll)^ 



we showed that the two threefold hollow sites, that are 
present on this basal surface, are energetically by far the 
most stable adsorption sites. Specifically, these are the 
fee site, located above a palladium atom of the third 
layer, and the hep site, directly above an atom of the sec- 
ond layer. Concerning O incorporation into the Pd(lll) 
surface, we will focus initially on the high-symmetry in- 
terstitial sites directly below the topmost Pd(llll) layer, 
postponing the discussion of deeper sites to the end of 
this paper. There are three such interstitial sites, namely 
an octahedral (octa) and a tetrahedral (tetra-I) site di- 
rectly below the fee and hep on-surface sites, respectively, 
as well as a second tetrahedral (tetra-II) site located be- 
low a Pd atom of the first substrate layer. The octa site 
is coordinated to six palladium atoms, while each of the 
tetrahedral sites is fourfold coordinated. The distinction 
between the two tetrahedral sites is due to the presence 
of the surface; in the bulk, they would be equivalent. 

The existence of three sub-surface sites together with 
two energetically favorable threefold hollow sites for on- 
surface adsorption accounts for a vast host of possible 
structural combinations that have to be considered. This 
is exemplified in Fig. ^ which shows all on-surface/sub- 
surface site combinations resulting for a total coverage 
of Otot — 0.75 ML within our set of (2 x 2) geometries, 
and one of the three O atoms located below the surface. 
Corresponding structures (not shown) for a total cover- 
age of 8tot = 0.50 ML (2 O atoms, 1 of which below the 
surface) and of Otot = 1.00 ML (4 O atoms, 1 of which 
below the surface) were tested as well. Very often two 
symmetry inequivalent geometries exist for the same on- 
surface/sub-surface coverage and site occupation, as is 
also illustrated in Fig. ^ We considered both cases, but 
found that in some situations one of them is energetically 
very unfavorable: For both the fcc/octa and hcp/tetra- 
I site combinations one of the possible structures con- 
tains a sub-surface oxygen atom located directly below 
an on-surface oxygen atom at a close distance. Due to 
the ensuing strong repulsion between the electronegative 
oxygens (cf. Ref. 0), these geometries always turn out 
significantly less stable than all other possibilities. Elec- 
trostatic repulsion between the oxygen atoms seems also 
to play a significant role for the energetic preference of 
all other structures. For each pair of symmetry inequiva- 
lent site combinations, the structure in which the oxygen 
atoms are located further away from each other (within 
the constraints of a (2 x 2) arrangement, cf. Fig. |2Jl is 
found to be the more stable one. The coordination of- 
fered by the sub-surface site is another decisive factor, 
since in most cases the tetrahedral sites are favored over 
the octahedral one. 

Within the set of structures with mixed on- 
surface/sub-surface site occupation the fcc/tetra-I geom- 
etry resulted then as the most stable one at all considered 
coverages (Otot = 0.50 ML, 0.75 ML or 1.00 ML). While 
initially pure on-surface adsorption is preferred, already 
at a total coverage of Otot = 0.75 ML the binding en- 
ergy of the mixed (2x2)- (2 0fcc + Ototra-i)/Pd(lll) 
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fcc/tetra-I (a) fcc/tetra-I (b) fcc/octa 



fcc/tetra-II (a) fcc/tetra-II (b) 




hcp/octa (a) hcp/octa (b) 



hcp/tetra-I 



hcp/tetra-II (a) hcp/tetra-II (b) 




FIG. 1: Top view of all possible on-surface/sub-surface site combinations for a total coverage Otot = 0.75ML in a (2 x 2) 
arrangement, and one of the three oxygen atoms below the surface. Inequivalent geometries with the same site occupation are 
denoted with (a) and (b), the (a)-geometries being the more stable ones. Big spheres correspond to palladium (white = surface 
layer, light grey = second layer), small spheres to oxygen (gray = on-surface, black = sub-surface). Sub-surface oxygen atoms 
in tetra-II sites are invisible in this plot and are indicated schematically by small white circles. 



TABLE L Binding energies, Et, (in eV/atom) for the on- 
surface adsorption in fee sites and pure sub-surface adsorp- 
tion into the three interstitial sites as a function of coverage. 
Additionally shown are the binding energies of the mixed on- 
surface/sub-surface occupation structure with O in fcc/tetra- 
I sites. Positive binding energies indicate that adsorption is 
exothermic. 
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structure exceeds the one of the (2 x 2) — 3 Ofcc/Pd(lll) 
geometry containing the same number of oxygen atoms. 
Since fee is the most stable adsorption site for on-surface 
adsorption^, this means that at the latest at 0.75 ML to- 
tal oxygen coverage, structures with oxygen incorporated 
into the sub-surface region become the favored ones. Fig- 
ure |21 gives a graphical summary of these results (cf. Ta- 
ble [}, where the dashed line marks additionally a Gibbs 
construction, i.e. the common tangent. In the coverage 
range spanned by the resulting shaded region (where the 
common tangent yields the lowest energy) we expect co- 
existence of the (fee/-) and the (fcc/tetra-I) structures. 
This coexistence may result in different domains of these 
two phases, or it may result in a mixed phase built from 
the elements of both structures. Obviously the coverage 
range of coexistence will be enlarged with temperature, 
and as said above, surface oxides will start to form as 
well. A quantitative evaluation of this phase transition 
and the corresponding actual composition of the surface 
requires an analysis using, e.g., Monte Carlo simulations 
(see e.g. Refs. il4Bl5il . which is, however, outside the 



FIG. 2: Average binding energy Et as a function of the to- 
tal oxygen coverage, with on-surface oxygen in fee and sub- 
surface oxygen in tetra-I site. The highest Eb is found for the 
pure on-surface chemisorption phase at Otot = 0.25 ML, but 
for 0tot > 0.50 ML the on-/sub-surface structures become 
more stable. A coexistence of the on-surface and on-/sub- 
surface structures is expected in the grey shaded region (s. 
text). The marks denote the actually calculated binding en- 
ergies listed in Table |I] The lines are drawn to guide the eye 
and are schematically extrapolated to zero. The dashed line 
is the Gibbs construction, i.e. the common tangent. 

scope of the present paper. 

Apart from coexistence, the determination of the 
threshold coverage Oc for oxygen penetration from the 
results summarized in Table |l] and Fig. [21 has also to be 
seen in light of the limited set of (2 x 2) surface unit cell 
structures considered. A quantitative determination of 
Oc would require a much larger set of structures (peri- 
odic and disordered) . For the present study, we therefore 
only conclude that Qc is roughly at 0.5 ML total cov- 
erage. Interestingly, the coverage above which bulk ox- 
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ide formation becomes thermodynamically stable is also 
~ l/2MLfi& and surface oxide formation will even set in 
earlie»ii. We believe that this points at an important 
bottleneck function of the initial O incorporation in the 
oxidation sequence of this surface: Sub-surface O can ap- 
parently not exist as a stable phase in the coverage range 
below which oxide formation would already set in, which 
in turn implies that the role of sub-surface oxygen is that 
of a metastable precursor. At not too low temperatures 
the phase transition to oxidic structures will then pro- 
ceed nearly instantaneously after the first incorporation 
of oxygen below the surface. 

A rationalization for this low stability of sub-surface 
oxygen can be found by analyzing the geometric struc- 
ture of configurations containing sub-surface O in more 
detail. The atomic geometry of the two most stable struc- 
tures for Qtot = 0.75 ML and Qtot = 1.00 ML is displayed 
in Figs. 12^^ andOt. Due to the presence of sub-surface 
oxygen lowering the symmetry, these geometries exhibit 
a much more complex relaxation pattern compared to 
structures with pure on-surface chemisorption, cf. Ref. 
0. However, a number of similarities in the relaxation 
pattern of different atoms in the two structures can be 
discerned. The three first-layer palladium atoms in the 
vicinity of the sub-surface oxygen move sideward away 
from it and are lifted up, increasing the volume which 
the sub-surface O atom has at its disposal. The fourth 
palladium atom is pushed away and downward, as the 
on-surface oxygen atoms in the immediate neighborhood 
of the sub-surface oxygen try to increase their distance 
to the latter, moving towards an on-top position. The 
slight change in the position of the other on-surface oxy- 
gen atom seems then rather as a consequence of this dis- 
placement of the surrounding palladium atoms. 

The vertical relaxations shown in Fig. Ofor both the 
(2 X 2) - (2 0fcc + Otctra-i)/Pd(lll) and the (2x2)- 
(3 0fcc + Otctra-i)/Pd(lll) system are calculated with 
respect to the center of mass of each layer. In addi- 
tion, there is a buckling of Azi — 0.28 A for the first 
and Az2 = 0.36 A for the second layer, respectively, 
and identical for both geometries. The most promi- 
nent feature is, however, a pronounced first layer ex- 
pansion, which is slightly smaller for the lower coverage 
case, i.e (2x2)- (2 0tcc + Otctra-i)/Pd(lll), reflecting 
the reduced on-surface oxygen coverage. Interestingly, 
the bondlength between the sub-surface O atom and its 
Pd neighbors remains almost constant for both relaxed 
structures and takes values between 1.9 A and 2 A. The 
corresponding numbers found for on-surface O atoms in 
our preceding study ^ were in the same range. In this 
light, we interpret the described relaxation patterns of 
the mixed on-surface/sub-surface structures largely as a 
consequence of a local expansion of the crystal lattice to 
accommodate the sub-surface oxygen at an optimal 0-Pd 
bondlength. 

Comparing these results for the initial penetration of 
oxygen into the Pd(lll) surface with corresponding pre- 
vious DFT studies at the basal surface of the left and 



a) 0.50 Ofcc/ 0.25 Otetra-I 




FIG. 3: Top view and side view of the most stable structure 
at (a) Gtot = 0.75 ML, i.e (2x2)- (2 Ofcc 4- Otctra-i) and (b) 
Otot = 1.00 ML, i.e. (2 X 2) - (3 Ofcc + Otctra-i). The small 
spheres represent the oxygen atoms, while the palladiums are 
depicted by the large spheres. Atoms, which are lying in the 
same plane and are equivalent under the present mirror sym- 
metry have the same color. The arrows (not drawn to scale) 
indicate the direction of the displacement of the substrate 
atoms from their ideal (bulk-like) positions. 



right neighbors of Pd in the periodic tableiSiiSiSfli^i, we 

are able to identify some clear trends"'^^. Concerning 
the critical coverage at which incorporation commences, 
occupation of sub-surface sites starts at Ru(OOOl) and 
Rh(lll) approximately at a full monolayer adsorbed at 
the surface, while at Ag(lll) this happens already at 
about just a quarter of a monolayer. We have showni^ 
that this clear trend is consistent with the bulk prop- 
erties of these late Ad transition metal elements, i.e 
it scales roughly with their bulk cohesive energies or 
bulk moduli. The already mentioned association to 
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0.00 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 

Oxygen coverage (ML) 

FIG. 4: Dependence of the calculated work function (left) 
and the dipole moment (right) on the total oxygen coverage 
for the most stable geometry at each considered coverage. 
At Gtot = 0.75 ML and Gtot = 1.00 ML both the fee and 
the fcc/tetra-I structures are shown. Neither the coexistence 
phase, nor surface oxides are considered here. 



coverages at which oxide formation has been observed 
also experimentalljiSiSSiSiS^iS^, further suggests that sub- 
surface oxygen acts as a metastable precursor to the for- 
mation of oxides at all of these metals. The sub-surface 
oxygen incorporation goes in all cases hand in hand with 
an appreciable deformation of the substrate lattice. In- 
corporation of 1 ML oxygen in a tetra-I site causes e.g. 
relaxations of -h53% on Ru(OOOl), -^48% on Rh(lll), 
+46% on Pd(lll) and -f42% on Ag(lll)i2ii2i2a*2i. Such 
a distortion is very costly for a comparably hard material 
like Ru, while it becomes less so for the more noble metals 
towards the right of the transition metal series, i.e. to- 
wards silver. This is also consistent with the pronounced 
tendency towards sub-surface island formation computed 
for Ru(0001)i'5 and Rh(lll)2a, which we find to be much 
weaker for Pd and Ag2i. The stability of oxygen in a 
sub-surface interstitial site is therefore mainly character- 
ized by minimizing the repulsive interactions with other 
nearby O atoms, in competition with a reduction of the 
energy cost for the substrate lattice deformation achieved 
by locating the O atoms close to each other. Sharing the 
high cost for the lattice expansion leads in Ru(OOOl) and 
Rh(lll) to clustering of the oxygen atoms in the sub- 
surface region. In Pd(lll) and Ag(lll) the smaller cost 
reduces this tendency, but leads still to a lower stability 
of sub-surface O compared to on-surface adsorption. 



B. Electronic structure 

Turning towards the electronic structure, we first inves- 
tigate the coverage dependence of the work function and 
compare the cases of pure on-surface chemisorption to the 
mixed on-/sub-surface site occupation geometries. The 



change in work function upon increasing oxygen coverage 
relative to the clean surface (^^^i^^" = 5.25 eV) is shown in 
Fig. ^ As expected for an electronegative adsorbate, ac- 
comodating oxygen on the surface strongly increases the 
work function, while the incorporation of oxygen into the 
sub-surface region clearly lowers the work function com- 
pared to the pure on-surface adsorption structures. This 
difference remains almost constant, ~ 0.2 eV for all mixed 
fcc/tetra-I structures, and reflects the competition aris- 
ing between sub-surface and on-surface oxygen atoms for 
the bonding charge of the surface Pd atoms. This renders 
the on-surface O species slightly less negatively charged 
(as can be seen by comparing the initial state shift of 
the O Is core levels), with respect to a situation with no 
sub-surface oxygen involved. The concomitant depolar- 
ization can be quantified via the dipole moment /i, which 
is related to the adsorbate-induced work function change 
by the Helmholtz equation 

Ai= (l/127r)AA$/e. (1) 

Here the surface dipole moment is measured in Debye, A 
is the area per (1x1) surface unit cell (in A^), the induced 
work function change, A<I'(8) = cl>(Q) — cjjcioan^ jg given in 
eV, and the coverage 8 in ML. As apparent from Fig. 01 
a strong and inward pointing dipole moment is obtained 
for all considered coverages, consistent with the strong 
polarization of the density at the surface oxygen atom(s). 
Concerning the pure on-surface adsorption the observed 
reduction of the dipole moment with coverage is then a 
consequence of the growing repulsive interactions within 
a more and more densely packed overlayer, which induces 
a slight depolarization of the O atomaS. The decrease of 
/i for the mixed on-surface/sub-surface geometries, can 
be understood along similar lines. Now, the polarization 
goes in the other direction and yields a dipole moment 
contribution of the opposite sign, reducing the overall fi 
of the surface. 

To gain further insight into the bonding mechanism un- 
derlying the superior stability of the Ofcc/Otctra-i struc- 
tures, it is helpful to also analyze the difference electron 
density. This quantity is obtained by subtracting from 
the electron density of the 0/Pd(lll) system both the 
electron density of the clean Pd(lll) surface and that of 
the isolated oxygen layers, each time keeping the atoms 
fixed at the same positions. Figure |S1 (left) exemplifies 
this for the (2x2)- (2 0fcc + Otetra-i)/Pd(lll) geome- 
try. A strong depletion of the Mxz,yz states of the sur- 
face palladium atom(s) bonded to both Ofcc and Otctia-i 
is discernible, while the electron density at the oxygen 
atoms is polarized. Comparing with equivalent plots for 
other on-surface/sub-surface structures, it appears that 
the linear Ofcc ^ Pd — Otctra-i coordination furthers the 
strong hybridization of the Pd4:dxz,yz orbital with the 
Ofcc and the Otetra-i species. Simultaneously, the Pd 
atom in between effectively screens the electrostatic in- 
teractions between the two O atoms, resulting overall in 
the high stability of this structural arrangement. 

Interestingly, a similar linear arrangement 0-TM-O 
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FIG. 5: Plots of the difference electron density for the (2 x 2) ~ 
(20fcc + Otctra-i)/Pd(lll) structure visuaUzing the bonding 
within (left) and the coupling of the O-Pd-0 trilayer to the 
Pd(lll) substrate (right). The contour plot depicts the [211] 
plane perpendicular to the (111) surface. Regions of electron 
accumulation and depletion have positive and negative signs, 
respectively. Contour lines are drawn at 0.06 A~"^ spacings. 
The nomenclature Oi and O2 refers to the two on-surface 
oxygen atoms marked as "1" and "2" in Fig. El 

(TM — Ru, Rh, Ag) was also observed in the struc- 
tures with mixed on-/sub-surface site occupation that 
were reported as most favorable at the basal surface of 
the neighboring elements of Pd in the periodic system of 
elementsi2i22iSi. For both the fcc-structured Rh and Ag, 
the most stable mixed on-/sub-surface site occupation ge- 
ometry also consists of a fcc/tetra-I site combination22i2i. 
For ruthcnium.^^, a hep material, initially a similarly hn- 
ear hcp/octa structure seems to be the most stable one. 
A barrierless displacement along the [112] direction leads 
to an energy gain of about 0.2 eV and results then also 
in a fcc/tetra-I arrangement of the oxygen atoms, with 
the first Ru layer lying in a stacking fault position. This 
comparison indicates that the mechanisms playing a role 
towards stabilizing one mixed site structure over another 
are quite similar for all these late Ad transition metals. 

It has been suggested that an improved coupling to the 
underlying substrate is ultimately the reason for the pref- 
erence of the fcc/tetra-I combination over the hcp/octa 
oneiS. The difference electron density plot for palladium 
(Fig. O also indicates the presence of such a strong cou- 
pling between the 0-Pd-O "trilayer" and the underlying 
palladium atoms. However, this coupling is even better 
visualized by another type of difference electron density. 
Here the electron density of the trilayer and the under- 
lying Pd(lll) substrate are subtracted from the over- 
all density of the full structure, as exemplified for the 
(2 X 2) - (2 Ofcc + Otctra-i)/Pd(lll) structure in Fig. 
(right). The apparent coupling, which seems to be medi- 
ated predominantly by the hybridization of the VdAd^-^ 
orbital with the O 2s orbitals of the sub-surface oxygen 
atom, can be also quantified by 

771 1 / niall 77iO-Pd— O I7'sub^ (r,\ 

^coupling — 2^-^ ^ ~ ) . (Z) 

Here i?**" is the energy of the whole system, £;0-Pd-o jg 



the energy of the trilayer and E^^^ is the energy of the un- 
derlying Pd(lll) substrate, as described above. The divi- 
sion by two accounts for the two surfaces present in our 
symmetric slab geometry. With this we determine the 
coupling of the trilayer to the substrate as 127meV/A^ 
in the (2 x 2) - (2 Ofcc + Otctra-i)/Pd(lll) structure and 
as 142meV/A2 in the (2 x 2)- (3 Ofcc + Otctra-i)/Pd(lll) 
structure. In the other on-surface/sub-surface structural 
combinations, the coupling is much weaker, especially for 
the alternative linear arrangement with hcp/octa site oc- 
cupations. In the (2 X 2) - (3 Ohcp + Oocta)/Pd(lll) the 
coupling is e.g. only 56 meV/A^, i.e. less than half of the 
value in the fcc/tetra-I structure at the same coverage. 
After maximum distance between O atoms and a linear 
on-surface/sub-surface arrangement, this confirms the 
coupling to the underlying substrate as another infiuen- 
tial quantity determining the stability of these structures. 
The "coupling difference density plot" in Fig. [S] (right) 
reveals furthermore that most of this coupling is the re- 
sult of the interactions between the sub-surface O atoms 
and the second layer Pd atoms, while the first layer Pd 
atoms (i.e. the ones within the trilayer) contribute only 
little. Coupling to the underlying substrate has also been 
identified as playing an important role in the stabiliza- 
tion of the surface oxide structures on both the Pd(lll)^ 
and the Pd(lOO)^^ surfaces. For these two, as well as for 
the surface oxide identified on the Ag(lll) surfaceSLS^, 
oxygen atoms are similarly present at the interface to 
the underlying substrate. In the case of Pd(lOO) the cou- 
pling has been calculated as lOOmeV/A^^ i.e. of the 
same order of magnitude as now found for the most sta- 
ble fcc/tetra-I geometries. 



C. Bulk dissolved oxygen 

So far we have investigated the presence of sub-surface 
oxygen between the first and the second layer of the 
Pd(lll) surface. It is conceivable that oxygen can also be 
incorporated between deeper lying layers, i.e layers which 
are further away from the surface. The incorporation of 
oxygen into the bulk material, on which we will focus in 
the following, is the limiting case for such a penetration, 
and will furthermore allow us to estimate the amount of 
dissolved O atoms in a palladium crystal. To address this 
question the binding energy of an oxygen atom in an octa- 
hedral interstitial site has been calculated. Considering 
the much smaller volume in the tetrahedral interstitial 
sites and the above described strong lateral expansion of 
the crystal lattice to achieve optimum sub-surface 0-Pd 
bond lengths of the order of 2 A, we expect the tetra- 
hedral sites to be significantly less stable (and thus less 
relevant for the total uptake) than the octahedral ones. 
Concentrating therefore on the octahedral sites, we em- 
ployed a large (4x4x6) palladium bulk unit cell with 97 
atoms and allowed relaxation of the nearest Pd neighbor 
shell around the interstitial site. Due to the still finite 
size of the supercell long range elastic interactions of the 
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metallic lattice are not properly accounted for in the re- 
sulting binding energy of Ei, — — 0.70eV. In a similar 
study concerning the oxygen concentration in a ruthe- 
nium crystals the long range elastic interactions were 
estimated to increase the calculated binding energy by no 
more than 0.5 eV. The calculated bulk modulus of ruthe- 
nium {Bq = 320 GPa) is nearly twice as large as that 
of palladium {Bq = 163 GPa). This leads us to assume 
that a consideration of the long range elastic interactions 
will probably not improve the calculated binding energy 
by more than 0.30 eV, leading to a conservative upper 
limit for the binding energy of bulk dissolved oxygen of 
Et « — 0.4eV. With respect to gas phase O2, occupation 
of interstitial sites is therefore endothermic. Still, a fi- 
nite concentration can nevertheless be expected at finite 
temperatures purely on entropic grounds. Together with 
the vast number of available interstitial sites, a consider- 
able amount of oxygen might therefore still be deposited 
into a macroscopic sample, which we estimate with the 
following thermodynamic reasoning: 

In thermodynamic equilibrium the number of oxygen 
atoms, A^O; that are dissolved in the bulk, depends on 
the total number of available interstitial sites, N, in the 
sample, as well as on the temperature T and pressure 
PO2 in the surrounding. By minimizing the Gibbs free 
energy one obtains for the concentration of oxygen in the 

Pd buiya 



^ [-Bi.+AMo(T.po2)]/fcBT 

N 



(3) 



where is the Boltzmann constant, Eb the binding en- 
ergy with respect to half of the total energy of molecular 
oxygen (positive for exothermicity) and the gas phase 
conditions are summarized in the chemical potential, 
A/xo(T,po2)i which is also referred to 5O2. The latter 
quantity is related to gas phase temperature and pressure 
via ideal gas lawa^S 



AMo(r,_po2) = ^ 



(4) 



where the temperature dependence of /io2 iT,p'^) includes 
the contributions from rotations and vibrations of the 
molecules, as well as the ideal gas entropy at p° = 1 atm. 
It is listed in thermodynamic tables^ or can be easily 
obtained from first-principles calculations^. 

Choosing a pressure of — 10~^^atm, character- 
istic for ultra-high vacuum (UHV) conditions, and us- 
ing the above derived value for the binding energy, the 
Nq/N concentration is obtained as 6.2 x 10~^* and 
1.2 X 10-1-* for T = 300K and T = 800 K, respectively. 
The number of octahedral interstitial sites in fee Pd is 
6.5 X 10^^ sites/cm"^, which translates the above concen- 
trations into 0.4 X 10^ O-atoms/cm"^ for a temperature of 
300 K and to 0.8 x 10^ O-atoms/cm^ for 800 K. Yet, most 
often, surface science experiments quantify uptakes in 
units of monolayers, which makes it desirable to convert 
the above numbers also into ML for direct comparison. 
Unfortunately, this requires to use a specified volume of 



the crystal. For convenience we take this to be Icm"^, 
which is of the same order of magnitude as a typical ex- 
perimental sample. Then a coverage of 1 ML corresponds 
to 1.5 X lO^^atoms/cm^ and the above stated concen- 
trations translate into the equivalents of 10~^°ML and 
10~^ML at room temperature and T — 800 K, respec- 
tively. We conclude that the amount of bulk-dissolved 
oxygen in a macroscopic Pd crystal in equilibrium with 
its UHV surrounding is negligible. 

Of course, during adsorption experiments, the crys- 
tal surface is exposed to much higher pressures than the 
UHV pressure considered above. On the other hand, 
it is unlikely that during the finite exposure time, the 
whole sample attains already the saturation concentra- 
tion corresponding to full equilibrium with the surround- 
ing higher gas pressure. If we therefore consider am- 
bient pressure in our thermodynamic estimate, it will 
correspond to a conservative upper limit: The effec- 
tive pressure during exposure will not exceed it, and the 
real uptake in the finite exposure time must be smaller 
than the equilibrium value. However, even when using 
PO2 — 1 a-tni as the other extreme for the estimate in eq. 
©, we only obtain about lO^^'ML (10" atoms/cm^) 
and 0.6 ML (10^^ atoms/cm"^) dissolved oxygen at room 
temperature and T — 800 K, respectively. We therefore 
always end up with uptake numbers that are significantly 
smaller than the fewi to few hundred^ ML of oxygen 
reported in experimental studies. Even considering the 
uncertainty due to the employed gradient-corrected DFT 
functional, the discrepancy between the experimentally 
reported oxygen uptakes and our calculated ones seems 
simply too large. We are therefore forced to conclude 
that the only place to store such large amounts of oxy- 
gen atoms in a Pd crystal is not in the regular intersti- 
tial sites considered here, but in bulk defects like grain 
boundaries and dislocations, or in oxidic structures pre- 
sumably forming in or at the surface of the crystal. This 
conclusion is therewith equivalent to the one of a preced- 
ing study on the dissolution in Ru, where even smaller 
uptake numbers were computediS. Considering that Pd 
is a much softer and larger lattice constant material with 
a correspondingly reduced lattice deformation cost upon 
O incorporation, the larger amount of oxygen dissolved 
in palladium is not astonishing. Unfortunately, corre- 
sponding studies have not yet been done for Rh and Ag, 
but one would expect the number of oxygen atom solved 
in Rh to be intermediate between the ones found for Ru 
and Pd, whereas the palladium value will be exceeded by 
Ag. 

Returning finally to the initial motivation of studying 
bulk dissolution as opposite limit to the immediate sub- 
surface incorporation studied in sections HI A and B, 
we notice that the endothermic bulk binding energy of 
Eh ~ — 0.4eV is comparable to the binding energies cal- 
culated for pure sub-surface oxygen adsorption, cf. Ta- 
ble I, i.e. for incorporation below the surface without 
adsorbed oxygen present on the surface (-/tetra-I). Both 
are therewith significantly less stable than adsorption in 
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the mixed on-surface/sub-surface structures, leading us 
to conclude that O accommodation occurs preferentially 
in the immediate selvedge of the crystal. The ensuing 
strong local enrichment with oxygen entails then natu- 
rally the formation of ordered oxidic structures at the 
surface. 



IV. SUMMARY 

We presented a density-functional theory study ad- 
dressing the initial oxygen penetration into the Pd(lll) 
surface. The results show that sub-surface oxygen is ini- 
tially always less stable than on-surface oxygen adsorp- 
tion. However, with increasing coverage this preference 
vanishes, and oxygen incorporation becomes more favor- 
able already above coverages of '--^ 1/2 ML. This threshold 
coverage is similar to the coverage for which bulk oxide 
formation becomes thermodynamically stable-'^^. and for 
which surface oxide formation was reported experimen- 
tally. We interpret this as reflecting a bottleneck func- 
tion of the initial O penetration in the oxide formation 
process, with an almost immediate phase transition to 
oxidic structures after the first incorporation of oxygen 



below the surface. In the resulting structures with oxy- 
gen in on-surface and sub-surface interstitial sites, we 
always observe a clear preference for linear 0-Pd-O ar- 
rangements, with occupation of fee and tetrahedral sites. 
Similar trilayer structures had already been reported in 
preceding studies on the O incorporation into the basal 
surfaces of the other late Ad transition metals Ru, Rh 
and Ag. Analyzing the geometric and electronic struc- 
ture of the fcc/tetra-I trilayer, we identify a good hy- 
bridization and screening within the trilayer, as well as 
an improved coupling to the underlying substrate as key 
factors towards the stabilization of this geometry. Com- 
paring with bulk dissolved oxygen, our calculations indi- 
cate a strong enrichment of incorporated oxygen in the 
immediate selvedge of the crystal. Via the formation of 
trilayer structures, this suggests that the enrichment will 
then naturally entail the formation of oxidic structures 
at the surface. 
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